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ABSTRACT 

The lens mass determined from the photometrically obtained Einstein time scale suf- 
fers from large uncertainty due to the lens parameter degeneracy. The uncertainty 
can be substantially reduced if the mass is determined from the lens proper motion 
obtained from astrometric measurements of the source image centroid shifts, S9 C , by 
using high precision interferometers from space-based platform such as the Space Inter- 
ferometry Mission (SIM) and ground-based interferometers soon available on several 
8 - 10 m class telescopes. However, for the complete resolution of the lens parameter 
degeneracy it is required to determine the lens parallax by measuring the parallax- 
induced deviations in the centroid shifts trajectory, AS0 c . 

In this paper, we investigate the detectabilities of S9 C and AS6 C by determining 
the distributions of the maximum centroid shifts and the average maximum deviations 
expected for different types of Galactic microlensing events caused by various masses. 
From this investigation, we find that as long as source stars are bright enough for 
astrometric observations it is expected that S9 C for most events caused by lenses with 
masses greater than 0.1 M & regardless of the event types can be easily detected from 
observations by using not only the SIM (with a detection threshold Sdth ~ 3 //as) 
but also the ground-based interferometers (with <50 t h ~ 30 /tas). However, detection 
of AS6 C from ground-based observations will be difficult for nearly all Galactic bulge 
self-lensing events, and will be restricted only for small fractions of disk-bulge and 
halo-LMC events, for which the deviations are relatively large. From observations by 
using the SIM, on the other hand, detecting Ad6 c will be possible for majority of disk- 
bulge and halo-LMC events and even for some fraction of bulge self-lensing events. For 
the complete resolution of the lens parameter degeneracy, therefore, SIM observations 
(or equivalent) will be essential. 
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1 INTRODUCTION 



The biggest difficulty in revealing the nature of Galactic 
dark matter by using microlensing is that one cannot deter- 
mine the individual lens masses because the Einstein time 
scale, which is the only quantity related to the lens mass ob- 
tained from the photometrically measured event light curve, 
results from the combination of other lens parameters by 



AGM D i(D 0S — D i) 

D^s J ' W 

where M is the lens mass, v is the lens-source transverse 
speed, and D a i and D os represent respectively the distances 
to the lens and the source from the observer. As a method to 
resolve the lens parameter degeneracy for general microlens- 
ing events, it was proposed to measure the shifts of the 
source star image centroid caused by lensing from follow- 
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up astrometric observations for events detected photometri- 
cally from the ground by using high precision interferometers 
from space-based platform such as the Space Inter ferometry 
Mission (SIM) and ground-based interferometers soon avail- 
able on several 8 - 10 m class telescopes such as the Keck 
and the Very Large Telescope (Miyamoto & Yoshii 1995; H0 
g, Novikov, & Polarev 1995; Walker 1995; Dominik & Sahu 
2000). 

When the lens-source separation is normalized by the 
angular Einstein ring radius #e, the centroid shift vector 
with respect to the position of the un-lensed source for an 
event caused by a point-mass lens is represented by 



t-t 

t-E 



(2) 



(3) 



where u is the lens-source separation normalized by 6e, u x 
and u y represent the components of u parallel and normal to 
the lens proper motion vector fi, to is the time of the closest 
lens-source approach, and u m i n is the lens-source separation 
(also normalized by Oe) at that moment (i.e. the impact 
parameter). The angular Einstein ring radius is related to 
the physical parameters of the lens by 



vtv _ ( 4:GM D os ~D ot 
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(4) 



When u — \/2, the amount of the centroid shift becomes 
maximum with 



71' 



(5) 



During the event, the trajectory of the centroid shifts traces 
an ellipse (astrometric ellipse). Once the trajectory is con- 
structed from astrometric observations of the event, one 
can determine 6e because the size of the astrometric ellipse 
(semi-major axis) is directly proportional to 6>e- Then one 
can determine the lens proper motion by fi = fe/iE com- 
bined with the photometrically determined value of tE- With 
the determined value of /i, the uncertainty in the determined 
lens mass can be substantially reduced because fi does not 
depend on v. However, the lens parameter degeneracy is not 
completely broken even with the determination of n, because 
it still results from the combination of the lens mass and the 
location. 

The lens parameter degeneracy can be completely bro- 
ken, and thus the lens mass can be uniquely determined if 
one determines the lens parallax by astrometrically measur- 
ing the deviations in the observed centroid shift trajectory, 
A<50 c , caused by the change of the observer's location during 
the event due to the orbital motion of the Earth around the 
Sun: parallax effect.n However, A.56 c , in general, is much 



* The lens mass can also be determined by astrometrically mon- 
itoring nearby very high proper motion stars working as a lens 
instead of monitoring lensed source stars because the distance to 
the lens is known. This method was suggested by Miralda-Escude 
(1996) and Paczyhski (1998) and the specific candidate stars from 
the SIM monitoring were selected by Salim & Gould (2000). We 
note, however, that our analysis is focused on followup observa- 
tions for events detected by the photometric lensing surveys from 
the ground. 
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Figure 1. The parallax effect on the astrometric behavior of an 
example halo-LMC event. The upper left panel shows the source 
star trajectories with respect to the lens (L) with (solid curve) 
and without (dotted curve) the parallax effect. The dashed cir- 
cle centered at L represents the Einstein ring. In the upper right 
panel, the trajectories of the centroid shifts S0 C resulting from the 
individual source trajectories in the upper left panel are marked 
by the same line types. In the lower panel, we present the tra- 
jectory of the deviation vector A50 c , which corresponds to the 
vector difference between the solid and dotted curve in the upper 
right panel. Both the trajectories of S0 C and A<50 c are presented 
in the ecliptic coordinates. For the details about the lens system, 
see the text. 



smaller than SB c , and thus measurements of the parallax- 
induced deviations require higher precision observations. 
The possibility of astrometric parallax measurements was 
investigated by Paczyhski (1998) and further analysis con- 
sidering the actual performance of the specific instruments 
was done by Boden, Shao, & Van Buren (1998). However, 
their analyses were based on several example events. 

In this paper, we determine the distributions of the 
maximum centroid shifts, /(<5f? Cimax ), an d the average maxi- 
mum deviations, /((A<5# Cjmax )), for different types of Galac- 
tic microlensing events caused by various masses. From the 
analysis of these distributions, we statistically investigate 
the detectabilities of S0 C and A<50 c expected from ground- 
and space-based observations. 



2 PARALLAX-INDUCED DEVIATIONS IN 
CENTROID SHIFT TRAJECTORIES 

Due to the orbital motion of the Earth around the Sun, the 
lens location in the ecliptic coordinates changes by 

Aipx = nsin(A Q -A), (6) 
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Atpp = — II sin f3 cos(A0 — A), 



(7) 



where II = 1 AU / D a i represents the lens parallax, (A, 0) 
is the ecliptic coordinates towards the direction where the 
event is occurred, and \q represents the ecliptic latitude 
of the Sun during the event. Then the change of the source 
location (normalized by 9&) with respect to the lens becomes 



Au x 
Au v 



A(fi\ A<pg . 
— — cos <p H — — sm < 



( Pl, 



f D os -Doi \ 

V Dos J 



(8) 



( 



Alfa \ f D os - D l 

sm <p H — —— cos <p J 



Oe 



where <j) represents the angle between fi and the ecliptic 
plane. The factor (D os — D i)/D os is included because all 
lengths are projected onto the source plane, and the minus 
sign is included because the source position changes towards 
the direction opposite to that of the lens position change. 
Then the parallax-induced deviation vector of the observed 
centroid shift from the centroid shift expected without the 
Earth's orbital motion is obtained by 

asg c = se c (u + Au) - se c (u), (io) 

where Au = (Au x , Au y ). 

The parallax effect on the astrometric behavior of a 
microlensing event is demonstrated in Figure 1. The upper 
left panel shows the source star trajectories with respect to 
the lens (L) with (solid curve) and without (dotted curve) 
the parallax effect. In the upper right panel, the centroid 
shift trajectories resulting from the individual source tra- 
jectories in the upper left panel are marked by the same 
line types. In the lower panel, the trajectory of the devi- 
ation vector AS9 C is presented, which corresponds to the 
vector difference between the solid and dotted curves in the 
upper right panel. The event is assumed to be caused by 
a lens located in the Galactic halo with D D i — 10 kpc and 
M — 0.5 Mq for a source star located in the Large Magel- 
lanic Cloud (LMC) with D 03 = 50 kpc. We also assume that 
the lensing parameters are U m ln = 0.2 and £e = 100 days, 
and the closest lens-source approach occurs at the moment 
when the source is at the largest angular separation from 
the Sun (i.e. A© — A = 180°). Observations are assumed to 
be carried out during — 5£e < iobs < 5£e with respected 
to the time of maximum amplification. The ecliptic latitude 
of a star in LMC is f3 ~ sin _1 (-0.99). Note that both the 
trajectories of S0 C and AS8 C are presented in the ecliptic 
coordinates. From the figure, one finds that AS8 C is much 
(nearly an order) smaller than 89 c . 



3 DETECTION OF CENTROID SHIFTS 

In this section, we show that as long as source stars are 
bright enough for astrometric observations^ one can detect 
S8 C for most Galactic events caused by lenses with masses 



t Note that astrometric observations of microlensing events are 
restricted by source brightness. With the SIM having a detection 
limit of V ~ 20 mag, although most Galactic bulge events can 
be observed, a substantially fraction (nearly half) of the LMC 
events cannot be observed owing to the faintness of their source 
stars (Han & Chang 1999). 
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Figure 2. The distributions of the maximum astrometric cen- 
troid shifts expected from different types of Galactic events caused 
by various masses. Note that the distributions are arbitrarily nor- 
malized. 



greater than 0.1 Mq regardless of the event types from ob- 
servations by using not only the SIM but also the ground 
interferometers. 

To show this, we compute the distribution of the max- 
imum centroid shifts by 



/(<5#c,max) = / dD 0S p(D 



oo poo 



dD ip(D i)irr E 



< I I dv x dv y vf(v x ,v y ) 
'o Jo 

6 69c ,max 



(11) 



where te = D i9 E is the physical radius of the Einstein 
ring, p(D i) and p(D oa ) represent the matter density dis- 
tributions of the lenses and source stars along the line of 
sight, 8 represents the dirac delta function, v x and v v are 
the two components of v, and f{vx,v y ) is their distribution. 
For the physical and dynamical distributions of matter in 
the Galactic disk and bulge, we adopt the models of Han & 
Gould (1995). For the matter in the Galactic halo, we adopt 
an isothermal sphere model with a core radius of the form 
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Figure 3. The distributions of the average maximum devia- 
tions in the centroid shift trajectory for different types of Galactic 
events caused by various masses. Note that the distributions are 
arbitrarily normalized. 
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Figure 4. The distributions of the probability of detecting the 
parallax-induced deviations in the centroid shift trajectories for 
given values of detection threshold, S9 t h, as a function of (loga- 
rithmic) lens mass. 



^m=po4^4. (12) 

where r is the distance measured from the Galactic center, 
Ro — 8 kpc is the adopted Galacto-centric distance of the 
Sun, and po = 7.9 x 10 -3 Mq pc~ 3 is the halo mass density in 
the solar neighborhood. The adopted velocity dispersion and 
the core radius are a = 250/ km s _1 and r c = 2 kpc. The 
factors 7rr| and v in equation (11) are included to weight 
the event rate by the cross-section (i.e. Einstein ring area) 
and the transverse speed. 

In Figure 2, we present the computed distributions 
/(5#c,max) for different types of Galactic eventsj^j In each 

t We note that both the distributions /(<50 c ,max) and 
/({A56 c max )) are determined for all events with u m i n < 1.0, 
not considering the dependence of the detection efficiency on the 
event duration and the amplification. As a result, they are differ- 
ent from the distributions for events actually detected by the pre- 
vious and the current lensing surveys (Alcock et al. 1993; Aubourg 
et al. 1993; Udalski et al. 1993; Alard & Guibert 1995; Abe et al. 
1997). We determine the distributions in this way because the 
detection efficiency is greatly dependent on the adopted observa- 



panel, we mark the type of events according to the order 
of the lens and source locations. For example, 'disk-bulge 
events' refer to those caused by lenses located in the Galac- 
tic disk for source stars located in the bulge. From the ob- 
tained distributions, one finds that for a given lens mass 
bulge self-lensing (bulge-bulge) events are most likely to 
have the smallest value of <5# c ,max- One also finds that the 
expected most frequent value of <5# c ,max even for this type 
of events caused by small mass lenses with M ~ 0.1 Mq ex- 
ceeds ~ 30 /xas. Then, considering the detection threshold is 
60th ~ 3 ^as for the SIM (Unwin, Boden, & Shao 1997) and 
59th ~ 30 (j,as for the ground interferometers (Colavita et al. 
1998), detection of 56 c will be possible for most of events 
regardless of the event types from observations by using not 
only the SIM but also the ground interferometers. 



tional strategy, the instruments, and the data analysis technique, 
which are evolving rapidly, and thus the analysis based on the ob- 
servation condition and strategy of a specific lensing survey will 
have little meaning. 
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4 DETECTION OF PARALLAX-INDUCED 
DEVIATIONS 

In the previous section, we show that the amount of the max- 
imum centroid shifts for most Galactic microlensing events 
are large enough to be detected both from space and ground- 
based astrometric observations. However, as shown in the 
example of Figure 1, the amount of the parallax-induced 
deviations will be much smaller than the shifts caused by 
the lensing itself. Then a naturally arising question is that 
for what fraction of events one can determine lens parallaxes 
by detecting A.50 c . In this section, we answer this question. 

To determine the fraction of events with detectable 
A<5# c , wc compute the distribution of the average maximum 
deviations by 

poo pD os 

/((A50 c , ma x» = / dD os p{D os ) I dD ol p(D ol )nrl 
Jo Jo 

f 00 f°° (13) 
x / / dv x dv y vf(v x ,Vy) 

Jo Jo 

x S [<A5(9 c , max ) - (A4,m«)' (0 E ,Il;D o ,,D O3 ,M)] . 



event 


lens mass 


detection 


probability 


type 


(M ) 


SIM 


ground 


bulge-bulge 


0.1 


0.11 


0.00 




0.3 


0.27 


0.00 




0.5 


0.37 


0.00 




1.0 


0.53 


0.00 


disk-bulge 


0.1 


0.79 


0.26 




0.3 


0.86 


0.36 




0.5 


0.88 


0.41 




1.0 


0.91 


0.47 


halo-LMC 


0.1 


0.86 


0.10 




0.3 


0.93 


0.21 




0.5 


0.94 


0.28 




1.0 


0.96 


0.39 



Table 1. The detection probabilities of the parallax-induced de- 
viations in centroid shift trajectories which are expected from the 
astrometric observations by using the SIM and ground interfer- 
ometers. The assumed detection thresholds are 58^ = 3 fias for 
the SIM and S6 th = 30 ^tas for the ground interferometers. 



For the computation of (A<5# c , max ), we first compute A59 c 
of each event for given values of the lens parameters (M, 
D i, and D os ) and the parameters of the relative observer- 
lens-source positions (A0 and cj>) by using the equations in 
§ 2. We then obtain the maximum value which is expected 
from the assumed duration of observation. Since astromet- 
ric observations will be carried out for events detected from 
photometric monitoring, we assume observation duration to 
be —He < tobs < KMe- Since the lens proper motion vec- 
tor \i has no preferred orientation, the average value of the 
maximum centroid shifts is obtained by assuming that (f> is 
randomly distributed in the range < <j> < 2n. In addi- 
tion, since we are interested in all lensing events not in the 
events detectable by a specific lensing survey, the impact 
parameters of events have random distribution in the range 
< itmin < 1-0. The ecliptic latitude of a star towards the 
Baade's window is (3 ~ — 6°.5. For the physical and dynami- 
cal distributions of the lens and source, we assume the same 
distributions used for the computation of /(<5# Cjmax ). 

Figure 3 shows the obtained distributions /((A(5# Cjmax )) 
for various types of Galactic lensing events. In Figure 4, we 
also present the distributions of the probability of detect- 
ing A5# c for given values of detection threshold, 88^, as 
a function of (logarithmic) lens mass. From the figures, one 
finds the following facts. First, if observations are performed 
with the ground interferometers, it will be difficult to detect 
A.S8 C for nearly all bulge self-lensing events since the de- 
tection probability is less than 1% even for events caused 
by lenses with M = 1.0 Mq. In addition, detection will 
be limited for only small fractions of disk-bulge and halo- 
LMC events. For example, for these types of events caused 
by M = 0.3 M Q , the probabilities are 36% and 21%, re- 
spectively. Second, if events are observed by using the SIM, 
on the other hand, it is expected to detect the deviations 
for majority of disk-bulge and halo-LMC events and even 
for some fraction of bulge self-lensing events. For example, 
the detection probabilities for events with M = 0.3 Mq are 
86% and 93% for disk-bulge and halo-LMC events and 27% 
for bulge self-lensing events. In Table 1, we summarize the 



detection probabilities which are expected from the obser- 
vations by using the SIM and ground interferometers. 



5 SUMMARY AND CONCLUSION 

We compute the distributions of the maximum astromet- 
ric source star image centroid shifts, <5# c ,m ax , and the aver- 
age maximum deviations in the centroid shift trajectories, 
(A<5# c , max ), expected for different types of Galactic events 
caused by various masses. The findings from the analysis of 
these distributions are summarized as follows. 

(i) As long as source stars are bright enough for astromet- 
ric observations, one can detect SG C for most events caused 
by lenses with masses greater than 0.1 Mq regardless of the 
event types from observations by using not only the SIM but 
also the ground interferometers. 

(ii) If observations are performed with the ground inter- 
ferometers, detecting A.80 c will be possible for nearly none 
of bulge self-lensing events and for only a small fractions of 
disk-bulge and halo-LMC events. 

(iii) However, if events are observed by using the SIM, 
detection of A.S8 C will be possible for most disk-bulge and 
halo-LMC events and even for some fraction of bulge self- 
lensing events. 

Therefore, for the complete resolution of the lens degeneracy 
and accurate lens mass determination, SIM observations (or 
equivalent) will be essential. 
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